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ABSTRACT 


Studies have been conducted of the frequency stabilization 
of the output from a Sylvania 941 Carbon Dioxide laser by 
active comparison with a Stark-tuned NH.D line center in an 
absorption cell. A dithered Stark voltage gave a peak- 
detection correction signal fed back to a PZT mounted laser 
resonator mirror to control the cavity optical length. 
Operation of an integrator incorporated in the feedback 
loop to improve the system stability has been investigated. 
A method was demonstrated of using the integrator to locate 
the desired laser vibronic transition by scanning the laser 
through its spectral signature prior to initiation of servo 
operation. An unexpected discontinuity in the absorption 
lineshape ‘of a Stark-tunable cell was discovered which made 
the cell unsuitable for generation of the feedback signal. 
A new and more versatile Stark cell was designed for test 


and evaluation of cell plate configuration and fabrication. 
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I. INTRODUCTION 


The Sylvania Model 941 CO. Laser to be used in atmospheric 
extinction studies was determined to have poor beam quality 
and frequency stability of the laser output. Since the 
chance of oscillation on a specific vibrational-rotational 
line is dependent on whether a cavity mode is in resonance 
with a vibrational-rotational peak, the mechanism used 
in frequency tuning a co. laser involves varying the cavity 
mode spacing by changing the cavity length. As the laser is 
tuned in frequency by changing the mirror separation, 
competition effects occur and oscillation can switch from 
one rotational level to another. The fitting of a piezo- 
electric driver to one of the laser cavity mirrors has made 
possible the use of active feedback control to frequency 
stabilize the output. 

The use of active stabilization to minimize frequency 
drift has been approached in a number of ways. It requires 
a comparison of the laser output frequency with a suitable 
frequency standard, generating a correction signal to be 
fed back to the piezoelectric transducer to adjust the laser 
cavity optical length. Such a signal may be obtained by 
passing a fraction of the laser output intensity to a 
suitable detector through an absorption cell containing 
a suitable absorbing gas. The Stark effect permits tuning 


an absorption line of NH5D to oscillate about the desired 


laser vibronic transition frequency. The detector then 


ii 





Gives a Sign~dependent error signal. The ability to 
generate a discriminant voltage without impressing any 
modulation on the output and the voltage control of the 
frequency reference make the Stark cell stabilization 
technique attractive as a control system. 

This thesis describes the continuation of studies of 
the operation of Stark cell modulation, and the technical 


requirements for stabilization of a carbon dioxide laser. 





Pie  CHEORETICAL 


A. CARBON DIOXIDE LASER 
1. Lasing States and Spectra of Carbon Dioxide 

The quantum states and energy levels on which 
molecular lasers operate originate from three contributions: 
i) electronic energy, due to the electronic charge cloud 
surrounding the nuclei; ii) vibrational energy, due to the 
motion (vibrations) of the nuclei about their equilibrium 
positions; and iii) rotational energy, due to the rotational 
motion of the molecule about various axes in space. Each 
of these components is quantized. Each widely spaced 
electronic level is split into a series of more closely 
Spaced vibrational levels. This energy level scheme is 
shown in Figure 1. Molecular gas lasers exploit transitions 
between the vibrational-rotational levels of a molecule. 

In general, an N-particle system has 3N degrees of 
freedom. The motion of the center of gravity and the 
Orientation of the axis of a linear XY, molecule can be 
specified in terms of five degrees of freedom, leaving 
(3N - 5 = 4) four vibrational degrees of freedom. For a 


linear triatomic symmetric molecule such as CO there are 


5! 
three distinct normal modes of vibration. The four vibrations 
correspond to the symmetric stretch mode, the doubly degen- 
erate bending mode and the asymmetric stretching mode. These 
modes are shown in Figure 2. In a) the symmetric stretch 


mode, the oxygen atoms vibrate along the internuclear axis 


9 
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in a symmetric manner. In b) the bending mode, the atoms 
vibrate symmetrically in a direction perpendicular to the 
internuclear axis. The bending mode has a two-fold degeneracy 
which arises from the equivalence of vibration in either of 
two orthogonal planes whose intersection lies along the 
molecular axis. Inc) the asymmetric mode, the atoms vibrate 
asymmetrically along the internuclear axis. Since carbon 
dioxide is symmetrical in its equilibrium state it does not 
have a dipole moment. The symmetrical stretching vibration 
does not create a net electric dipole moment and therefore 

is not active in the infrared. The asymmetric normal 
Vibrations produce a dipole moment and are I.R. active. 

The stretching vibrations have higher frequencies than the 
bending vibrations since the restoring force is greater in 
stretching a molecule than in bending. 


The CO. molecule vibrational quantum state can be 


2 


described by a set of quantum numbers written as 


Q 
itt v5 v3 | 


where v is the degree of excitation of each of the three 
individual normal modes. The degenerate bending vibration 


mode for co. has an associated angular momentum which 


depends upon the excitation and relative phase of the 
orthogonal vibrations. This angular momentum is denoted 


by the quantum number £ appearing above v The superscript 


2° 
Q takes the value g = Vos V5 -2,..-, l, O, with each 


Giving rise to a sublevel. 





Assuming independence of the three normal modes, and 
the use of normal coordinates in describing the molecular 
vibrations, allows the separation of the Schrodinger equation 
into as many Separate equations as there are normal modes. 
The eigenvalues for each of the separate equations are 


given by 


Ea Fa) nev, (1) 


a: 


where Vs is the quantum number associated with the jth mode, 
and the total vibrational energy associated with an 


Iv vy, v3 | State is approximately 


= (v, + 45) ixeAs, ae MA Ge 4) hev5 + (v3 + 43) hey; oe) 
where VY is the frequency in wave numoer units associated 
with the normal modes of vibration. 

The carbon dioxide vibrational energy levels of 
principal interest in laser action with their modes and 


frequencies are 


ie} [00° | Asymmetric Stretch, yy, 23493 aa 
2) [10°o | Symmetric Stretch eS Sr om 1 
3) [02°o] Bending, 1225.5) ene 
4) [o1*o] Bending B73) em 
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The two vibrational levels [10° | and [02°o] belonging to 
different vibrational modes have nearly the same energy 
and are in Fermi resonance. When the interaction between 
them is included as a perturbation, the wave functions are 
mixed and the resulting pair of mixed levels are repelled, 
so that the level separation is greater than expected. 
The lower vibrational energy level diagram for the 

CO» molecule in conjunction with the No molecule is shown 
in Figure 3. The dashed lines indicate radiative transitions 
that are allowed by electric dipole considerations. The 2, 
and tT notation represents the electronic-vibrational 
angular momentum of the molecule (4=0,11=1), giving the 
value of the vector sum of the electronic orbital angular 
momentum and vibrational angular momentum about the inter- 
nuclear axis. The + superscript implies that the vibronic 
eigenfunction is unchanged when reflected at any plane 
containing all the nuclei. the subscripts indicate whether 
the state is u, antisymmetric (odd symmetry) or g, symmetric 
(even symmetry) with respect to vibrational motion. (Ret. 1] 

| For every vibrational state there exists a distinct 
set of rotational levels with slightly different spacing 
for different vibrational levels. The CO, molecule can 
rotate about any axis normal to the internuclear axis. 
From quantum theory of rotational spectra and angular momentum 
the energy level due to rotational motion is given by 


Ep (J) = BI(JI + 1) = eG eae) (3) 


ens. eee 
8n~cI 


3.3 
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maere J is the angular momentum quantum number (J = 1, 2,...), 
I is the moment of inertia, and B is called the rotational 
constant. The selection rule for J for dipole transitions 
Ot co. allows only transitions between rotational levels 
BeeewhichAJ = 11. 

Laser action for carbon dioxide occurs on transitions 
between two vibrational levels of the same electronic state 
Or more precisely, transitions between rotational levels of 
mie two vibrational states (vibration-rotational transitions). 
These transitions are classified into "branches" by the 
change in J. Transitions between two vibrational levels v 
and yv' with account taken for the rotational splitting are 
shown in Figure 4. The transition, in the absence ot 
rotational energy, would consist of a single line centered 
at Vo! called the Q branch where 4J = OQ. Including rotational 
energy the transition actually consists of two groups of lines: 
the P branch corresponding to a change in rotational quantum 
number of AJ = +1 and R branch corresponding to a change in 
rotational quantum number of AJ = -l. 


Solving for Y in the equation 
QE = hy (4) 
and defining ve as the band center frequency or 


Voz E, (0) - By (0)] (5) 


cs 





the transition frequencies can be written as 


Ys ie Bay P(J + 1) Branch 
Vis = Ve Q Branch (6) 
yey + (*Byn)g R(J - 1) Branch 


where J is the rotational quantum number of the lower level. 
The rotational levels of a linear molecule are positive or 


negative depending on whether the sign of the total eigenfunction 
Sos = Tgugte Coe (7) 


remains unchanged or changed upon an inversion. 
The infrared selection rules for the vibrational- 


rotational transitions relevant to co. are [Ref. 2] 


> 
pen) 
ll 
© 


~-l; Gafeg; usu (8) 


AJ = 0, _1l; +++; Sexe a (J = 0)** (J = 0) 


where s and a designate the symmetric and antisymmetric 

Getal rotational eigenfunction J and the symbols «—~ and 

af» represent allowed and not allowed respectively. In all 
symmetric excited vibrational levels (x7) the even rotational 
levels are positive and the odd are negative with respect to 
inversion along the molecular axis while the converse is true 


for the antisymmetric levels (eee hor the fir, A... douply 
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Segenerate vibrational levels, for each J there is a positive 
and negative level with alternating order (+-, -+, +-, -+,... 
or ~-+, t-, -+, +-,...). 

Because co. 1s exactly symmetric and is without 
Bm@etear spin the QO branch or AJ = 0 transition is strictly 
forbidden by symmetry restrictions (7+—~ x , 43 #0). J to 
g 1 transitions within a Single vibrational energy level 
are also forbidden. Since the two oxygen atoms are spinless, 
symmetrically located indistinquishable particles and the 
CO, wave function remains unchanged, the system obeys Bose- 
Einstein statistics. Therefore all co. g states have 
rotational levels with even J values only, while all u states 
have only odd J values. The absence of alternating rotational 
levels causes the spacing between adjacent lines to be twice 
that expected on the basis of the size of the molecular 
rotational constant B. This spacing is 4B instead of 2B. 
Experimental evidence shows that only the even J values 
appear in co. P- and R-branch spectra. The “=e Cy, system 
can oscillate on approximately one hundred transitions 
Spanning the wavelength range from 9.1 to 10.7yum. They 
occur on the P(J) and R(J) branches of both the 00°1+—> 10°0 
and 00°1*~ 02°O rotational-vibrational bands. A detailed 
laser transition diagram for the lobe - 10°0 and 00°1 = 02°0 
bands including some rotational levels is shown in Figure 5 
and some of the P- and R-branch transitions for the 00°1 - 
10°0 vibration-rotation levels with their corresponding 


transition wavelengths and frequencies are given in Table I. 


[Ref. 3] 
a, 
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TRANSITION 


P(14) 
B16) 
P (18) 
P (20) 
P (22) 
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P (28) 


TRANSITION 


Bi4) 
R (16) 
R(18) 
R(20) 
R (22) 
R (24) 
R (26) 


R (28) 


TABLE I 


oO 
00-1 - 10°O TRANSITION 


P BRANCH 


WAVELENGTH ( m) 


HO 75326 
Oe 5516 
LOzD/23 
JAC oe EY 
10.6118 
10.6324 
10.6534 


10.6748 


R_ BRANCH 


WAVELENGTH ( m) 


10.2860 
LOe2S55 
10.2605 
10.2470 
OR 2335 
10.2200 
LOS 2Z0N > 


Oro OU 


IES, 


FREQUENCY (cm 
949. 
947. 
945. 
944.18 
942. 
940. 
938. 


O36: 
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78 
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2. Excitation Mechanism 

Electrical pumping ina gas laser is achieved by 
allowing a current to pass through the gas producing ions 
and free electrons. The average electron energy and the 
ability of an electron to excite a neutral molecule by 
collision is much greater than the ion energy. The equili- 
brium condition established among the electrons can be 
described by an effective electron temperature bs Assuming 
a Maxwellian energy distribution, and that for each collision 
some fraction § of the electron energy is lost, the effective 
electron temperature is given by [Ref. 4] 

earrs ee (EL) (9) 

where E is the applied electric field, k is Boltzmann's 
constant, e is the elementary charge on an electron, and g 


is the electron mean free path. Since the mean free path is 


,=- =, (20) 
is Cea 2 
with fino the collision cross section and P the pressure, 


T. can be given as a function of the ratio a P. The ratio 
of applied electric field to pressure is used as a parameter 
in describing the discharge condition. 

The primary electrical pumping mechanisms for a 
Sas are: i) excitation by electron impact or a collision 


of the first king represented by the equation 


20 





A + (e + KE) —> A* +e (11) 


where * represents the excited state of A, and KE the 
translational kinetic energy converted to excitation 
energy. ii) Excitation by collision between atoms or 
molecules of different species (resonant energy transfer) 


characterized by the equation 
A+ B¥ —»> A* + B + KE (12) 


where A* and B* are the excited states. This type of 
excitation process is called a collision of the second kind. 
Selective excitation of polyatomic molecules through 
vibrational energy transfer from a diatomic molecule has 
been used effectively in the production of population 
inversion on vibrational~-rotational transitions of the 
polyatomic molecule. For the N, - co. system shown in 
Figure 3, the [00°1] level of carbon dioxide is chiefly 
excited by a resonant energy transfer via collisions of the 
second kind between co. molecules in the ground state and 
N. molecules in the excited state. As a result of these 
resonant collisions the populations of certain sets of 
levels are very strongly coupled so that they act essentially 
as a Single level. The vibrationally excited N, molecules 
in the ground electronic state N,* (v = 1) are produced 


primarily by indirect electron impact excitation. The 


pump rate R,, of the upper state is 


BAL 





ve) 
tl 
Put 


p q <vo> (13) 


where 


<vo> { vo (£) £(E) dE (14) 

N is the electron density, V is the electron velocity, 
f(E) is the electron energy distribution and O(E) is electron 
excitation cross section. [Ref. 5] The cross section for 
electron impact excitation of N. in the ground state is 


very high making substantial population densities of N_* (v =1) 


2 
achievable. Figure 6 shows the total electron excitation 
cross section of the N. levels as a function of electron 
energy with the cross section reaching a maximum value of 
approximately 3 X 107 +© em? at an electron energy of 2.3ev. 

[Re£. 6] The vibrationally excited nitrogen molecule in 

the ground electronic state cannot decay through electric 
dipole radiation back to the ground state (v = 1 to v = Q) 

since N, is a diatomic homonuclear molecule with no permanent 
dipole moment. The energy difference E between the N.* (v = 1) 
and co, [00°1]’ energy levels is about 18cm +. ‘This is 
negligible in comparison with the average thermal energy 

KT at room temperature of 210em ~. The near resonance 

molecule - molecule collision leading to energy transfer 


can be written as 


O 2 O = 
N.* (v=1) + co. [00 o| --- N. (v=0) + co. [00 1] -18em (15) 


ac 





This resonance transfer produces a mixed state in which the 
combined populations of N.* (v = 1) and CO. [00° 1] are 
essentially in equilibrium: this effect increases the 
effective lifetime of co, [o0°1] . The transfer cross 
section is largely due to the negligible difference between 
the energy levels. The level [10°0] (or [02°] ) is more 
than 900em * from N,* (v = 1) and the transfer cross section 
is small compared to that for the reaction in (15). 

To a lesser degree excitation can occur by direct 


Swectron impact on the ground state CO, molecule [00° o] 


2 
and by resonant exchange of energy bdetween CO molecules 

in excited levels and co. in the ground state. Electron 
impact on the ground state CO, molecules results in excitation 
of the CO, molecule to the [00°21] Keve mcreniece ethic trans lt2on 
co a [oo° 1] is optically connected to the ground state 


.. [00°0} and is given by 


+ O 15 2 
The cross section for vibrational excitation of co. by 


electron impact is shown in Figure /7. [Ret. 7] 

A substantial amount of CO (10% of co., pressure) 
can be produced as a result of dissociation of — ina 
discharge. The electron excitation cross sections for CO 
are reported in Reference 8 to be similar in form to that 
of N. but shifted to lower energies with the peak value for 


2 
; -16 
the total cross section being approximately 8 X 10 cm 


2 3 





at 1.8 ev. The difference between the first CO vibrational 
energy level and the upper laser level of co. is given as 
igen *. These final two mechanisms account for the 
possibility of producing a population inversion in a pure 
CO5 laser. 

Normally the co. oscillator functions with a plasma 
that consists of a CO., = N, ~ He mixture. The helium 
component plays two important roles. The first is related 
to the increased rate of relaxation of the [0170] level. 

By reducing the bending mode population which in turn reflects 
itself in a reduction of the lower laser level populations 
the gain is increased. Helium also reduces the kinetic 


and rotational temperature of the CO, in the plasma. This 


2 
will also increase the gain, since the same total upper- 


state population is now distributed among fewer levels. 
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3- «Lineshape Function 
The function describing the distribution of emitted 
intensity versus the frequency, Y, is called the lineshape 


function, g(Y) and is normalized according to 
6} 
(gv) ay = 1 (17) 
=03 


The frequency spectrum and g(v) are the same for 
purposes of spontaneous emission, stimulated emission and 
absorption. Line broadening mechanisms may be divided into 
two classifications: i) homogeneous and ii) inhomogeneous. 

Homogeneous broadening involves the spectral spread 
of the whole system caused by the broadening of the spectrum 
of each indistinguishable molecule. It is due most often 
to the finite collisional lifetime of excited molecules 
and the consequent interruption of the radiative process. 
The result of frequent collisions is the presence of many 
truncated radiative or absorptive processes with a line- 
width of radiation greater than that of an individual 
uninterrupted process. 

The normalized Lorentzian lineshape function which 


is characteristic of homogeneous broadening is [Ref. 9] 


AVY. 


= Sa a 18 
2T[(v-vy)* + (AV. /2)" J ial 


gly) = 


Here Y, is the center frequency and 4Y is the separation between 


the two frequencies at which the Lorentzian is down to half 
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its peak value, or the full width at half maximun. 
The linewidth Avy, is related to the finite lifetime T 


of the emitting state by 





The Lorentzian linewidth is pressure dependent 
Since the frequency of collisions is proportional to the 
density of the gas. The time interval between collisions 
can be approximated by the ratio of the mean free path 


(10) to the mean thermal velocity VY 





= 8kT 4s 
v= ( TM ) (20) 
and 
L. 
aoe MKT) “¢ 
a (21) 
Cc (167 )2 Pa 


where M is the mass of the molecule. Since T is inversely 
proportional to the pressure P, the linewidth AV is 
proportional to the pressure. 

Inhomogeneous broadening is due to the spread in 
the individual transition frequencies of the molecules and 
not the broadening of the lines of the individual molecules. 
The thermal motion of the gaseous molecules is the cause 
of Doppler broadening. The transition frequency Y is 


Doppler shifted due to the finite velocity of the molecule. 
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Meeng the Maxwell velocity distribution function of a gas, 
the normalized Doppler broadened lineshape function for 


inhomogeneous broadening is Gaussian and given in Reference 10 


ao 


where M is the molecular mass at equilibrium at temperature T, 
and V . is the resonant frequency. 


The linewidth obtained from (22) is 


AY, = 2V0 (23) 
and reexpressing g(v) in terms of AY, gives 
2 ~ -Vv : 
ee exp | 40 2) (V0) 2 2 | (24) 


ee AV, 


For a carbon dioxide laser Doppler broadening 
dominates at low pressure giving a GausSian line, and 
pressure (collisional) broadening dominates at high pressures 
giving a Lorentzian line. In general, the overall broad- 
ening is due to a combination of effects and mechanisms and 
is given by a convolution of the various processes. The 
overall lineshape (go) of two broadening mechanisms, one 
homogeneous (g,) and the other inhomogeneous Sa) may be 


written as 


1 


g, W) = f a, (v-v'a,(v') av’ (25) 
= 00 
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The convolution of a single Lorentzian line with 


a Single Gaussian in (25) gives a Voigt profile. [Ref. 11] 


If both Doppler and collision widths are appreciable, the 


total line width can be approximated by [Ref. 12] 


AY 


Figure 
of Lorentzian, 


the lineshapes 


5 


= J av? vie 


] (26) 
Doppler collision 


8 (a) and 8 (b) shows a general comparison 
Gaussian, and Voigt lineshapes. In 8 (b) 


were normalized so that g(o) = 1 and in 


8 (a) the curves were drawn with equal linewidths. The 


Voigt profile was calculated numerically on a Hewlett- 


Packard 9810A calculator using (18) and (24) in (25). The 


program used was a slightly modified version of the one 


found in Reference 13. 
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4. Laser Modes and Spacing 
The general expression for a laser resonator cavity 


resonant frequency is given by 


Cc 


y eee SE 


(27) 


where 
1 1 L i "5 
OS Ga Ses f(a = Ry. (ees =| 


5 and R, are the radii of curvature of the two laser mirrors 
Seperaced by a distance L, 7 is the index of refraction of 
the medium, 2 and m are the transverse mode numbers with 
values of O, 1, 2,..., and q is the longitudinal mode number 
which is equal to the number of half wavelengths in the 
distance L. 


From (27) the separation between laser oscillation 


frequencies in the longitudinal (or Axial) mode for Aq=1 is 
AVq=Vq+ 1) - Vq = ©/2nL (28) 


If the gain linewidth of the medium is greater than or 
equal to the axial mode frequency spacing 4VYq, several 
longitudinal modes may oscillate at a time. 

For each value of q but different values of 2 +m 
there is a set of transverse mode frequencies falling between 


those which result from a change of the longitudinal mode 


ah 





index q. The frequency separation between adjacent trans- 


verse modes is 


ane lg aera : Wiese ie Dae Co 
In most applications of gas lasers, operation in the lowest 
order (2= m= 0) or fundamental transverse mode is desired. 
This can be accomplished by inserting into the laser cavity 
a circular aperture of size such that it passes the funda- 
mental mode with little loss but blocks much of higher order 
transverse modes which extend further from the optic axis. 
Typical gas lasers tend to operate in the fundamental 
transverse mode because of the aperture represented by the 
discharge tube. 

For the case of a resonator with mirror reflectivities 

oa and €,, an average distributed loss constant % and the 
average loss per pass equal to AL - mM (Cy e,)? the passive 


resonator linewidth is given by [Ref . 14] 


1 
my o[ x = TES Im. ea sao, ] . c (1 - R) (30) 
Cc 271 7] 2 Ti a) 1e 


Where in the approximation R = HE ay)” =i anadenx = 0. 

In general the passive cavity resonance frequency ae 
does not coincide with the molecular line center frequency 
Vo! resulting in the oscillation at a frequency near ue 
but shifted slightly toward Yo This effect is called 


frequency pulling and is due to a shift in the index of 
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refraction near molecular resonance. The atomic dispersion 


"pulls" the laser oscillation frequency Y from the passive 


resonator value toward that of the atomic resonance V ; 
fe) 


Giemitrequency of oscillation for a high Q resonator cavity 


is given by 


ze _ _ AVe 
Vevi- (Vo -¥o) es (31) 
where AV is the total molecular linewidth due to Doppler 
and pressure broadening. 
For operation in the lowest order transverse mode 
where & = m= 0 (27) becomes 


Vag = [gq a5 an cas = vo (32) 


Emesticution of (30) and (32) into (31) gives for the 


frequency of oscillation 





E C (q+£) c% (1-R) Vo ¢ (IER)] 1 
V= [q + £] 2m L | are 12 27 7 XS 


Equation (33) shows that the laser oscillation 


frequency may be tuned by changing the cavity length. This 


is the fundamental means of stabilization techniques using 


a piezoelectric driven mirror in the cavity. 


che: 


(33) 





9. Frequency Instabilities and Causes 


The condition for stable resonator operation is 


given by the condition 
L L 
OS (1 - “/R,) (1 - “/R,) = 1 (34) 


The change in laser frequency due to the change in 


cavity length is given by [Ref. 15] 
qy =-y —— (35) 


For a system where the cavity mirrors are fixed to a cavity 
structure with a coefficient of expansion, & , the 
change in cavity temperature OT can change the length of 


the cavity resonator and thus the laser frequency by 
Se eee. te OK AT (36) 


Therefore changes in ambient temperature and acoustic 
vibrations in the cavity which modulate the cavity length 
are sources of instability which can contribute to both 
long and short term drifts in the cavity resonance frequency. 
The oscillation frequency of a cO., laser can also 
change by varying functions in the operating parameters 
Such aS pressure and discharge current. Mocker in Reference 16 
reported a 5 - 3 MHz/torr downshift fe frequency due to an 


increase in pressure and a .5 - .9 MHZ/torr upshift in 
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frequency due to an increase in discharge current. He 
concluded that the pressure dependent frequency shift is 
caused by a change in the overall index of refraction of 

the gas discharge and the current dependent frequency 

shift is attributed to the change in the index of refraction 
caused by the change in the electron density of the plasma. 
Ina cO., gas discharge the refractive index 7) is determined 
by all constituents in the plasma such as the neutral 


molecules and electrons: (Ref. 17] 


2 
y=1-3 2, : + 27% (37) 
; Ma, Eo Ww ala 


Where 7} . and me are the electron density and mass 
respectively: 7 Se and x, are the density and polarizability 
of the neutral molecules respectively. A significant re- 
duction of the refractive index may arise from dissociation 


of CO, molecules into components with a smaller electric 


2 
Dermittivity. 

Therefore it is important to operate a cO., laser 
with a well regulated low ripple power supply and closely 
controlled gas pressure and temperature. 

The frequency deviation Ay from the resonant 


frequency Y is therefore determined by the fluctuation in 


the refractive index A” and cavity length AL ana iS given by 


Av  _ _4L , 47 (38) 
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For a sealed “* laser the pressure fluctuations are 
negligible and the temperature and flow rate of cooling 


water are the dominant causes of instability in the 


laser. 
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B. FREQUENCY TUNING MECHANISM 
By bringing a cavity mode into resonance with a 

vibrational-rotational peak a co. laser can be tuned to a 
vibrational-rotational line. The cavity mode spacing SV. 
of a Single mode a laser dictates the chance of 
oscillation on a specific rotational level and changes 

as the cavity mirror separation is varied. The optical 
length of the Model 941 laser cavity is 45 cm, short 
enough to ensure single axial-mode operation and give a 
value for Sv. = Yo - Vgtl = © 12 Or 353Mnhz. Ine cavity 
linewidth calculated from (30) with ® = 0 and using the 
reflectivities given by Sylvania for the cavity mirrors, 
va mee0 and €, = -/0 is 19MHz. Each vibrational- 
rotational transition is lineshape broadened with a line- 
width AY. From (23) the Doppler-broadened linewidth 

for tiie co. transitions is approximately 50MHz. For the 
pressures and mixtures used in single transverse mode 
(small bore diameter) ss lasers (15 - 20 Torr) the 
LorentZian pressure broadening AV, can exceed the Doppler 
broadening. The value listed by Sylvania for Av, is LOOMHz 
giving a total linewidth using (26) of 112MHz. The absolute 
frequency separation of the vibrational-rotational lines 
Sapp is between 52 and 57GHz depending on the J number 
@emtne transition. Figure 9, although not drawn to scale, 
shows three P{J) lines and a set of cavity modes showing 


the frequency spacing involved. The successive order of 


the transitions is dictated by the absolute cavity length; 
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they occur in an arbitrary fashion giving each laser its 
own spectral signature. 

AS a co. laser is tuned in frequency by changing the 
cavity length, competition effects occur and oscillation 
switches from one rotational line to another. In the 
interrotational level competition region, two rotational 
levels compete for the total population of the same and 
higher vibrational levels. Mocker [Ref. 18 | reported that 
the competition effects are due to intensity-dependent 
anomalous dispersion arising from index and gain saturation 
and that the switching from one rotational level to an 
adjacent level occurs over a small finite frequency interval. 

While the laser oscillates on a given vibrational- 
rotational transition the population density or the popu- 
lation inversion on that particular transition is reduced 
by saturation effects until the single pass gain in the 
medium equals the single pass loss of the optical cavity. 
If the vibrational level lifetime (typical values 1 - 10msec) 
is considerably longer than the rotational relaxation time 
me usec) a Boltzmann distribution will describe the 
population density of the various rotational levels at 
the same temperature. Therefore if the gain on the 
oscillating transition is the highest the gain on the 
non-oscillating vibrational-rotational transitions starting 
from other rotational levels in the vibrational level will 
be less than the single-pass loss of the optical cavity. 


This means that once the laser transition starts 
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Oscillating, the gain on the other transitions falls below 
that required for laser oscillation with the other rotational 
levels coupling (due to collisions) their stored energy 
into the lasing transition and only one transition 
Oscillates. This is the situation of very strong com- 
petition effects. As the cavity mode is tuned away from 
the center of the lasing transition a point is reached 
where another transition has equal gain and a competition 
will be initiated. The intensity dependent dispersion 
represents a frequency pulling term such that a cavity 
resonance 1s pulled toward the line center frequency. 
Transitions of approximately equal strength have their 
competition region halfway between line centers, while 
unequal strength transitions compete closer toward the 
weaker transition. Figure 10 illustrates the competition 
between vibrational-rotational transitions and how the 
laser can be tuned in frequency by varying the cavity mode 
separation as a result of changing the cavity length. [In 
Figure 10(a) the cavity mode at Y,, , is in resonance with 
the P(J,) tases transition causing it to oscillate and P(J>5) 
is between the cavity modes. In Figure 10(b) the cavity 
length is increased slightly and the modes at Voz and / yo 
are in resonance with the P(J,) and P(J.) transitions 
respectively and competition can occur. If the gains are 
mutually exclusive then only one line would survive at the 
expense of the other. If the gain of each line is sufficient 


to sustain oscillation simultaneously then both laser lines 
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would coexist. In Figure 10(c) the cavity mode at Ym3 iS 

in resonance with the P(J.) transition causing it to oscillate 
while the P(J,) transition is between the cavity modes. This 
frequency tuning mechanism makes it possible to use a piezo- 
electric driver to actively control and stabilize the 


meeaguency output of the laser. 
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3 CO, LASER P(J) TRANSITIONS AND CAVITY MODES 
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Cs FREQUENCY STABILIZATION TECHNIQUES 

Although many of the instabilities contributing to 
frequency fluctuations in the co. laser can be reduced by 
external means (ie. control of discharge current by 
regulated power supply and control of system temperature 
with the aid of a separate temperature controlled heat 
exchanger) to frequency stabilize the laser actively 
requires the ability to change either the index of refraction 
of the laser medium or resonator cavity length. The laser 
cavity length can be accurately controlled by mounting a 
piezoelectric crystal which expands and contracts proportion- 
ally to an applied voltage. An error or frequency correction 
Signal that is proportional to and indicates the direction 
of the required length correction can drive the piezoelectric 
crystal and frequency stabilize the laser. A common feature 
shared by a number of approaches involves the basic principle 
of establishing a frequency reference either internal or 
external to the laser to which the frequency is locked. 
By applying a modulation signal to the laser, causing its 
oscillating frequency to "dither about a reference point, 
and using phase sensitive detection techniques, an error 
signal can be fed back to servo the laser frequency to the 
fererence. 

Another stabilization technique which has some useful 
advantages over the methods referred to above uses a frequency 
reference that is external to the laser with the discriminant 


obtained not by dither modulating the laser frequency about 
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the reference, but by modulating the reference about the 
laser frequency. The important advantage of this technique 
is that it avoids impressing modulation (AM or FM) on the 
laser output. A gas cell containing a molecular absorber 
Which is tuned into resonance with the desired CcO., laser 
line by a suitable D.C. electric field applied to the 
absorber provides the frequency reference. By super- 
position of an A.C. electric field, the absorption line 
center can be modulated about the lasing frequency and 
phase detected, producing the required correction signal 
back to the piezoelectric drive on the laser. Another 
advantage of this technique is the fact that it is voltage- 
controllable which means the laser frequency may be con- 
trollably tuned across the transition line width by 


programming the D.C. bias voltage applied to the Stark 


cell. 
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D. STARK CELL STABILIZATION 
1. Stark Effect 

Stark effects are the changes in the spectrum of 
a system due to an applied electric field. The rotational 
spectrum of a molecule Which has a dipole moment is modified 
in an electric field, since the field exerts torques on the 
electric dipole moment changing its rotational motion. 

The effect of an electric field on the rotational 
energy levels of molecules is given by quantum-mechanical 
perturbation theory. The Stark effect in a gas occurs When 
an external electric field is applied to the gas. The effect 
is characterized by the splitting of an absorption peak 
into several peaks by removing the degeneracy of the 
rotational energy levels. The application of the field 
in the absorbing gas results in the splitting of a line 
iacOo 2) + 1 lines. 

The first order Stark effect which is linear with 
the applied electric field is produced by applying the elec- 
tric field to a symmetric top molecule (molecule with three 
or more fold symmetry - molecule with two principal moments 
of inertia equal). A symmetric top molecule's principal 
advantage is the existence of a permanent electric dipole 
moment. The Stark splitting of the energy levels in molecules 


possessing a permanent dipole moment 1s given by [Ref. 19] 


yr ME K 


hAy = \—_—_———- 
7 CJ (te) 


cos © (39) 
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mmere 4( 15 Che molecular dipole moment and E is the electric 
field strength, O is the angle between 4 and E, J is the 
total rotational quantum number, M is the magnetic quantum 
number or the component of J along a fixed direction, and 

K is the quantum number describing the interaction of the 
two principal moments of inertia. The selection rules for 


yipration-rotational transitions are 


OT 


I 
© 


aay l 
AK = 0 For Parallel Bands 
AK = +1 For Perpendicular Bands 


AM = O For Stark Field Parallel to Laser 
Electric Field 


AM = +1 For Stark Field Perpendicular to 
Laser Electric Field 


For asymmetric molecules or those which do not 
have a permanent electric dipole moment, the expression 
for the Stark splitting of the energy levels is more 
complicated and involves second order perturbation calcu- 
lations including the neighboring states J’ =J +1 ana 
a = J- 1. Tne energy splitting goes as the square of the 
external electric field. The second order quadratic Stark 
effect is usually considered negligible in Stark modulation. 
[Ref. 20] 

For slightly asymmetric top molecules an inter- 
mediate type of Stark effect results. These molecules 


can have pairs of nearly degenerate levels and a first 
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order Stark effect can result for small molecular 


asymmetry and large fields. 
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2. Selection of Gases 

Changes in absorption produced by an applied 
electric field on the states involved in transitions as 
given by the vibrational-rotational selection rules may 
result in a modulation of the laser intensity by a gas 
containing an appropriate molecule. In order to achieve 
an appreciable Stark effect and large modulation J should 
be small. 

The selection of a suitable molecule for modulation 
of an individual laser is rather straightforward. The 
criteria to achieve a high modulation index are: 1) it is 
necessary to use molecules with a large dipole moment; 

2) slightly asymmetric molecules may be used provided the 
Stark effect is essentially linear at the values of applied 
electric field used for modulation; 3) a linear Stark 
effect in either lower or upper state or both is suitable 
and the overall effect is determined by combination of the 
Stark effect on both states; 4) the absorbing Stark sensi- 
tive line and laser line should be within one or two 
Doppler widths of the absorbing gas. This requirement is 
not valid in the pressure broadening case where the higher 
pressures result in broader lines, poorer resolution and a 
smaller index of modulation with equivalent fields. 

Investigations by Landman, Marante and Early [Ref. 21] , 
Jensen and Tobin [Reé. 22| and most recently Martin, Corcoran 
and Smith [Ref. 23] have reported many gases capable of 


Stark modulation of a co. laser. Table II lists some of 
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the absorbing gases used for Stark modulation and their 


corresponding — 


laser lines. 


TABLE If 


GASES USED IN STARK MODULATION OF 00°1 - 10°0 TRANSITION 


Gas 


Laser Line Modulated 


Poy 2 ls) a (20), RV1G),) R( 138) ,.R(20) 
P (22) 

P (14) 

R (14) 

Pia pee (lowe ( 20)", RLS) 

PLO (12) (14) Pte P( 20), P22) 


P(14), P(20), R(12) 
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3. NHoD Rotational-Vibrational Spectra 


Of the gases in Table II the microwave and infrared 


spectra of NH.D have been studied extensively. Brewer et. al. 


2 
in Reference 24 reported a rapid Stark tuning rate and 
relatively high absorption coefficient for NH,D making it 
possible to obtain substantial modulation with relatively 

low electric field strength and short path length. Kelly, 
Paancke and Feld [Ref. 25] observed the transitions which 
could be tuned into resonance with three of the 10.6ym CO, 
laser lines, namely P(20), P(14) and R(12). The NH.D 
transitions were identified as belonging to the type 

(vo, J) = (0, 4)4—>(1, 5). The assigned transitions are 
shown in Figure ll with the corresponding co. laser lines 
enclosed in brackets. The lower and upper levels, the Stark 
shift ‘os required to bring each NH,D transition into resonance 
with its respective laser line, laser line wavelength (em +) 
and zero field splitting (em +) for the three laser lines 
are given in Table III from Reference 24. Of primary 
interest in this discussion is the transition involved in 


absorption of the P(20) CO, laser line in the 10.44m band. 


2 


The NH,D absorption line involved belongs to the 


transition 
(v, he he = (O., 49a) ora uae 205) 


This transition is nearly in coincidence with but at a 


slightly lower frequency than the 944.19cem + P(20) line 
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of the CO, laser. The P(20) laser frequency exceeds the 

zero field splitting of the (O_. 4  )ow—»(1 , 5.__) 

transition by 1720 + 60 MHz. A relatively large and linear 
Stark effect results from mixing of the two nearly degenerate 
rotation-inversion levels (O., 454) and (yee 454) tay alee 
lower state of the transition. The subscript s(a) refers 

to a symmetric (antisymmetric) inversion state with respect 
to simultaneous exchange of all pairs of identical nuclei. 


The upper level (lo, 395) has only a negligible second 


maces Stark effect [Ref. 26] . 
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RELEVANT ENERGY LEVELS OF NH, D 
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TABLE III 
IDENTIFICATION OF OBSERVED NH,D TRANSITIONS 
LOWER UPPER dbo LASER FREQ. ZERO FIELD 
LEVEL LEVEL (MHz) (cm7!) SPLITTING (cm!) 
Om; 40y, Vos 44 =2300=60 949 43 949 40 
0.9 Ag, lar 5g5 1720460 944.18 944.14 
0,9 44, 142 5p, +2800%100 970.55 970.64 
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4. Error Signal Analysis 


The modulated output intensity from a Stark cell 
due to absorption by a single Stark component, exhibiting 


a linear Stark effect is given by 
Do ol Gags -[08 | (40) 


Here Y is the absorption coefficient of the gas as a function 
of light intensity and applied electric field, I, is the 
input intensity and 2 is the cell length. The line center 
absorption coefficient where the line is collision-broadened 


in the high pressure Lorentzian region is given by [Ref. 27] 


2 
SS Tee. ae _ 32 
eS Slr (yg 8) (4 


Meme A,1s its wavelength, 7T_, is the radiative lifetime 
of the vibrational-rotational transition, N, and Ny are the 


upper- and lower-level population, respectively, 35 and oy 
are the corresponding statistical weights and dy, is the 
homogeneous linewidth as given in (19). The pressure 
dependence of AY is cancelled since N,, and 7 are also 
proportional to pressure, and A(A,) is a constant. 

The absorption coefficient at line center for an 


inhomogeneous line in the low-pressure Doppler broadened 


limit is given by (Ref. Pull 


p, L 
IL ro 2 Im 2 _ S2 
A( A.) = - aa TZ wv, ( = ) ( 5 a N) (42) 





Where Av, is the inhomogeneous linewidth given by (23) and 
is a constant. The pressure dependence of NI and Ny cause 
A( Ae) to vary linearly with pressure. 

Plant and Abrams in Reference 27 plotted the line 
center absorption coefficient vs Stark cell pressure for 
NH,D as shown in Figure 12. Using Voigt profile fits to 
the line shapes they were able to get a relation for the 
line center absorption coefficient which was an excellent 
fit to the experimental points. This relation is 

A( A.) = a Ag* eines! (N_- ale ) (l-erfa) (exp[ : ) « (43) 
2 8% To, AV (W)e 2 9g, 1 pla J 
L 
Where a = ( 4Y,/ Av, ) ( 2m 2) : and erf is the error function. 

From Figure 12 at pressures above about 5 torr one 
can approximate the Stark cell absorption lineshape with a 
Lorentzian and A(A,) = Dace? a constant. The mod- 
ulator bandwidth reaches a limit, since the Stark effect 
changes character as the frequency of the perturbing field 
becomes comparable to the linewidth. Claspy et. al. [ Ref. 28] 
estimated that the modulator bandwidth increased with 
pressure according to a pressure broadening coefficient. 
Therefore the broader line at higher pressures required a 
greater field variation for the same absorption change. 


The absorption coefficient in (40) can be given by 


Y = ON, a) ie g(V) (44) 
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Where & is a constant dependent on gas pressure and g(y ) 
is the lineshape function of the absorbing gas. Assuming 
a Lorentzian lineshape and letting «& equal the pressure 
broadening coefficient C, times the pressure P, (44) can 


be written as 


AY 


= A( No) Cr 45) 
i MAX 27r[ (y-v,)°- (4W/2F J | 


Mmaere AY is the linewidth of the Stark cell, and CViove ) 
is the frequency offset of the absorption line center from 


the laser frequency. Let (V-v,) be defined as € where 
eo = Cee SE (46) 


€, is the offset frequency dependent on the offset of the 
absorption line from the laser frequency at zero field and 

S is the Stark coefficient or Stark tuning rate which relates 
the line shift to the applied electric field. 


Substituting (46) into (45) gives 





T= ACA) CP = oa 
MA 


x ° AY /2) 7 
aw[(e. + SE)* + ( 4¥/2)*] 
Therefore, for a fixed pressure in the cell and a set laser 
line, 89 the absorption coefficient is only dependent on the 
applied Field E . Impression of an A.C. voltage with an 
I e 
Plectric field in the form E = AE coswt on the applied 


field E will give f a time dependence by oscillating the 
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absorption line shape about the line center frequency 
which will act as a dither mechanism about Yy . A schematic 
of the dither mechanism for a Stark cell modulator is shown 
in Figure 13. The remainder of the work described in this 
thesis will refer to the P(20) - NH,D combination. 
Nussmeier and Abrams [Ref. 26] give the frequency 
offset between the absorption line center of piss Je and the 


P (20) 10.6 wm line center as a function of applied voltage as 


E€ = - + SE = -2042 + 0.143|M|E (48) 
Where € is given in MHz, M is the 2 component of the 
rotational quantum number J and has integer values from 
-4 to 4, and E is the applied electric field in “/en. 
Setting € = 0 and M = +4 in (48) shows that an electric 
field of 3570 “/cm is needed to bring one component of 
NH, D absorption into resonance with the P(20) transition. 
The values reported in the literature for the Stark tuning 
rate, S, for the P(20) co. line and NH.D with M = +4 are 


2, 
consistent with Nussmeier and Abrams value of .572 MHz Saye 


56 





DITHER MECHANISM 


ee 
Lee SIGNAL 
=—=—4 DETECTOR 


OSCILLATOR 








ABSORPTION COEFFICIENT 


1S(E. + AD cos wt) 
} 


| 
ee ace ena 
i 
{ 
t 





AY 


— wwe Ss SS —_ ff 22°62 eu 


FREQUENCY 


FIGURE, > 


a 





eee ocark Cell Stabilization System 

A schematic of a laser stabilization system is 
shown in Figure 14 and is a somewhat simplified version of 
the one given in Reference 29. A small percentage (less 
than me/om= to avoid saturation effects) of the laser 
Output is directed into the Stark cell via a beam splitter 
and plane mirror. The basic Stark cell consists of two 
parallel metal electrodes in a vacuum tight container with 
an entrance and exit window. A precision high-voltage 
power supply biases one Stark electrode to bring an 


absorption line in NH.D into resonance with the P(20) 


Z 
line center. An A.C. voltage from an audio oscillator 
drives the second electrode at a set frequency. The cell 
contains NH,D formed by mixing one part NH. with ND, to a 
desired pressure in the cell. The exchange of hydrogen 
and deuterium is reported in Reference 30 to be rapid with 
an equilibrium concentration to 45% NH,D reached in less 
than one minute. When the laser is oscillating on the 
P(20) line the Stark effect will produce a sinusoidal 
modulation in the intensity through the cell which can be 
detected by a pyroelectric or HgCdTe detector. The signal 
is phase sensitively detected relative to the reference 
audio oscillator waveform ina lock-in amplifier. The 
detected phase is dependent on the relative frequencies of 
the absorption line center and the laser output. The 


magnitude and sign of the lock-in amplifier D.C. output 


Which is proportional to the phase difference between the 
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laser frequency and the Stark-shifted reference will produce 
a suitable error signal. After this error signal is 
appropriately phased, integrated and amplified it is fed 
back to the piezoelectric mounted laser cavity mirror to 
tune out the frequency difference and the frequency control 
feedback system is complete. Because the laser is capable 
of oscillating on a number of vibrational-rotational laser 
transitions dependent on the resonator cavity length, a 
linear ramp voltage applied to the piezoelectric drive will 
cause the laser to scan through its spectral signature 
profile until the P(20) transition is reached at which time 


the servo stabilization technique can be initiated. 
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Tit. EXPERIMENTAL INVESTIGATION 


Preliminary engineering work done in Reference 13 
studied the characteristics of the Sylvania oo, Model 941 
Laser pertinent to the design of a Stark cell stabilized 
system. The instabilities of the laser due to fluctuations 
in discharge current caused by an unregulated, unfiltered 
laser power supply and system temperature variations were 
measured and discussed. The piezoelectric parameter and 
thermal expansion coefficient also were determined for the 
laser. Computer analysis demonstrated the feasibility of 
producing a suitable error signal for phase sensitive 
detection with available equipment and the need for a servo 
loop integrator. A vacuum tight Stark cell was constructed 
and tested to 3000 volts with no evidence of breakdown. 

The parameters for the Sylvania 941 piezoelectric 
tunable and linearly polarized CO. laser used in the 
following discussion of experimental procedures and results 
are given in Appendix A. Appendix B discusses the optical 
design and specifications of the optical engineering oo 
laser spectrum analyzer used in the experiments. 

To prevent the instabilities caused by changes in dis- 
charge current a stable high voltage power supply (Sylvania 
941S) which in turn was run from a stable line regulator 
was used to power the laser. The use of a circulation pump 


minimized flow rate fluctuations and the system usually 


oul 





reached an equilibrium temperature condition after a few 
minutes. The current to the laser from the 941S power 
supply could not be controlled manually and was preset 


for optimum performance of the laser. 
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A. SERVO LOOP INTEGRATOR 

In order to provide optimum frequency stabilization 
of the co. laser the P(20) line should be maintained in 
coincidence with the center absorption frequency of a 
Stark shifted eo transition as accurately as possible. 
The two transitions in coincidence will produce a zero 
error signal. As the P(20) line drifts in frequency the 
magnitude and sign of error signal will be proportional to 
the change in frequency, and direction of the drift, 
respectively. If a proportional correction signal is 
Supplied directly to the piezoelectric driver the magnitude 
of the error signal will decrease to an equilibrium condition 
at a nonzero value for €. The equilibrium condition exists 
when the PS voltage producing a frequency difference € is 
equal to the PZ voltage supplied by the correction signal. 

Integration of the error Signal will produce a correction 
Signal that will drive the piezoelectric crystal in a manner 
proportional to the total time integral of the error signal. 
This integration will allow the error signal to return to 
zero at the required P4 voltage. The integrated signal 
will continue to increase until the frequency difference 
(error signal) is zero. At € = 0 the integrator maintains 
mime voltage at a constant D.c. value until a nonzero error 
Signal indicates the need for a change. 

The integrator shown in Figure 14 serves a dual purpose: 
i) integration of the error signal supplied by phase detection 


in the lock-in amplifier when the stabilization feedback 
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Semerol System 1S in Operation. ii) Integration of a 
selectable D.C. voltage to provide a linear ramp voltage 
to the laser plezoelectric driver in order to search for 
the required laser transition before the servo stabilization 
procedure can be initiated. 

The integrator described in Appendix C was designed 
to respond to both slow and rapid changes in frequency of 
the P(20) line from the Stark shifted NH,D absorption center 
frequency. The gain was designed to be equal to one to make 
the integrator output equal to the integral of the input. 

A diagram of the apparatus used to examine the 
characteristics of the integrator 1s shown in Figure 15. 

The input to the integrator was provided by a regulated 
O - 25 V D.c. supply which was set at a constant voltage 
A three way reversing polarity switch S: was used to provide 
positive or negative values of V. to the integrator or to 
isolate the integrator from the input voltage. The input 
voltage to the integrator was fed into channel A of a 
time base plotter and the output of the integrator into 
channel B. Figure 16 shows the integrator output as a 
function of a step input with V; = +2 volts. The graphical 
output indicates that integration of a 2 volt input for 
approximately 10 seconds gives an integration output of 
20 volts and slope of 1.95 V/sec. Thus the integrator 
output equals the integral of the input. Figures 17 a &b 
respectively show typical operation and response of the 


integrator to step inputs of +4V and +6V. Since the 
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MaeegratOr 15 part of the feedback control system that 

drives the piezoelectric crystal, differences between the 
analog integration of the input signal and the exact math- 
metical values are not critical and can be compensated for 

in the proportionality constant associated with amplification 
Memes COrrection signal. In other words, the integrator 

is concerned primarily with the manner in which the piezo- 
electric crystal is driven rather than the magnitude of the 
Correction signal. 

The integrator saturates at a voltage of +14 volts; 
therefore if the P(20) line consistently drifts in the same 
direction the integrator may saturate. By ganging a polarity 
reversing switch to the input and output of the integrator, 
Saturation can be controlled by reversing the Sign of the 
Signal in magnitude while maintaining sign of the proper 
correction signal to tune out the frequency difference. 

Minor drift in steady state voltage output inherent in 
analog integrations can be controlled by a voltage offset 
null on the integrator. This possible integrator-induced 
long term instability will cause the laser frequency to 
drift slightly but it can be corrected for in the feedback 
control system like any other instabilities in the system. 

The incorporation of an integrator should improve the 
ability of the feedback control system to frequency stabilize 
the laser and produce a system suitable for laboratory or 


field use. 
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B. LASER SPECTRAL SIGNATURE 

The laser is capable of oscillation on a number of 
transitions other than the P(20) as the resonator optical 
length varies. This capability adds to the versatility of 
the laser as well as satisfying a need to search for and 
find the desired transition before the feedback stabilization 
system can be initiated. If the cavity length of a CO. 
laser is slowly increased a certain sequence of lines will 
appear. After changing the cavity length approximately one 
half wavelength the sequence will repeat itself. The output 
profile obtained by changing the laser mirrors separation is 
called the laser signature. 

A schematic of the experiment to determine the spectral 
Signature of the laser is shown in Figure 18. The spectral 
Signature was produced by plotting the relative intensity 
of the laser as a function of applied voltage to the 
piezoelectric driver (P4% voltage). The P% voltage was 
supplied by a Kepco model ABC regulated D.C. supply with a 
range of 0 - 1500 volts. It was operated in the voltage- 
programmable mode with a positive ground. It was pro- 
grammed to operate at a voltage ratio of 100:1 output to 
input. Using the integrator described in Appendix C to 
supply a linear ramp voltage as the input to the programmable 
power supply, the laser could be scanned through a free 
spectral range. The scan rate was controlled by the magni- 
tude of the input voltage to the integrator since the slope 


of the integrator output voltage is directly proportional 
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to the input voltage. The direction of scan was controlled 
by a 3 position reversing polarity switch that changed the 
Sign of the slope as a result of changing the sign of the 
input voltage. The output of the integrator which was 
directly proportional to the P% voltage was supplied to the 
X axis of an x - y plotter. The relative intensity of the 
chopped laser output was obtained from a pyroelectric 
detector signal that was amplified, filtered and rectified. 
The relative intensity signal was then supplied to the 

y channel of the plotter. The frequency of the laser for 

a given PZ voltage could not be monitored, but the 
vibrational-rotational transition could be identified 


using an Optical Engineering CO. spectrum analyzer described 


Z 
in Appendix B. Typical results are shown in Figure 19 for 
an increasing PS voltage scan (0 - 1500 V D.c.) and Figure 20 
for a decreasing P% voltage scan (1500 - 0 V D.C.). 

The power meter in Figure 18 monitored the operation 
of the laser. It was used as a comparison with the results 
shown in Figures 19 & 20. Figures 21 & 22 give typical 
results for relative power output as a function of applied 
PZ voltage. 

In each set of plots the scan direction appears to 
have caused a hysteresis effect [Ref. 24). Differences in 
the traces appear to be caused by hysteresis effects and 
nonlinearity of the piezoelectric element. If the dipoles 


(whatever their nature: electronic, ionic or polyatomic) 


resist the change of the electric field the electric 
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Gisplacement (strain) will lag the electric field (stress). 
A lag in the displacement relative to the voltage change 
results in a hysteresis effect and the scan direction will 
shift the laser signature in the direction of scan, indepen- 
dent of the scan rate. From Figures 19 - 22 nine laser lines 
were obtained in the laser signature and the free svectral 
range of the laser is approximately 1.4 kv which agrees with 
the value found in Reference 13. 

Looking at the shape of the laser Signature it is 
apparent that none of the lines are allowed to lase over 
their full width before the neighboring line begins to 
oscillate. The mode competition or the fact that the 
oscillations at various transition frequencies compete for 
the available gain are quite prevalent ina co. laser and 
help to insure single wavelength operation. Throughout 
the experiment only one vibrational-rotational line was 
observed to lase at a time. The laser oscillated in the 
lowest order transverse mode Beg throughout the full range 
of the PS voltage. 

“The method of scanning through the laser signature can 
be used in the stabilization system shown in Figure 14 to 
moma the required P(20) transition prior to initiating the 


feedback control. 


al 
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C. STARK CELL OPERATION AND TEST 
1. Stark Cell Design 

The Stark cell described in Reference 13 was 
expected to operate satisfactorily in the stabilization 
system. The cell was designed to eliminate a persistent 
problem of breakdown in previous cells. All metal surfaces 
were isolated from the low pressure gas filled channel 
through which the laser beam would pass. The cell was made 
to accomodate one inch NaCl windows at each end and had a 
Single port used to fill the cell with the NH_D gas. 
Fabrication of the cell resulted in an uneven plate separa- 
tion of approximately 7mm with 1mm of plexiglass between 
each plate and the channel. The cell was tested to 3000 
volts with no breakdown and held a good vacuum. Figure 23 
shows the cell design with pertinent specifications. 

As discussed in Section II. D. 3 the rest frequency 

of the (0, 4 ieee) ———- (1, 5 IM = 51) transition of 


04 05 
NH.,D is shifted approximately 1700 MHz into coincidence 
feenethe 10.6 ym P(20) line by applying a Stark field of 
3570 VY/em. The voltage required to produce the Stark field 


is given by 
V=- eo ag (49) 


where V is the difference in potential between the plates, 
E is the electric field strength between the plates and d 


represents the path of integration taken between the plates. 
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STARK CELL AND SPECIFICATIONS 
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Although the simplest path between the plates is a 
perpendicular straight line the electrostatic force is 

path independent since each plate is an equipotential 
surface. Applying (49) to Figure 24 for a straight perpen- 
dicular path from the lower plate (L) to the upper plate (vu) 


yields. 


V =-{ 5 -d a-{ 5 cos 180° d = os: af = E [d-(atb)|+E[Latb] 


Where Ea is the electric field with no dielectric, dis the 
separation between plates, a and b are the dielectric 
thickness at the plates and E is the electric field inside 


the dielectric. Since E can be related to Es by the relation 
ol ea 9S x (51) 


where K is called the dielectric constant, (50) can be 


written as 


x atb 
v=E [a - (ath) +e ( SS) (52) 
Substituting in the values for the Stark cell shown in 
Figure 21 and letting Eo = 3570 V/om. The approximate Stark 
cell voltage to bring one component of the center absorption 


frequency of NH.D into resonance with the P(20) line is 


2 
V = 2030 volts. Here d = 7mm, a = b = lmm and K = 2.7 - 3.2 


for plexiglass (methyl methacrylate). Although this is a 


TA 


(50) 





relatively high applied D.C. voltage, it is within the 


breakdown limits specified for the cell. 
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2. Stark Cell Absorption Lineshape 


Experimentation to determine the absorption lineshape 
and modulator performance as a function of pressure began by 
measuring the relative intensity out of the Stark cell as a 
Function of pressure with no electric field across the plates. 


The attenuation of monochromatic radiation follows Beer's law 


d == 
“x (vy) = - Aly) Iv) (53) 
from which 
-— AX 
I. = I, 2 (54) 


where 4 1S a linear extinction coefficient which is the sum 
of the coefficients of total absorption and total nonforward 
Beateerang. The relative value of 4 depends strongly on the 
density and molecular composition of the gas. Absorption is 
of greatest importance in the infrared (particularly in the 
ee) 24 wm wavelength). Rayleigh scattering because of its 
strong dependence on wavelength should be unimportant for 
wavelengths longer than the visible range and no Mie scattering 
should be present. Since the density is proportional to the 
pressure of the gas is dependent on the pressure. 

The experimental arrangement is schematically 
represented in Figure 25. The laser was tuned and maintained 
at the center of the P(20) transition by keeping the power 


meter at its maximum value and monitoring the spectrum 
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analyzer for the P(20) line. The cell was evacuated to 

27x ‘ea torr and used as a zero pressure reference. The 
Stark cell output intensity was measured by the pyroelectric 
detector. The y axis of the x - y plotter sampled the 
relative intensity of the chopped Stark cell output via the 
detector, amplifier and filtered rectifier. The x-axis of 
the plotter was manually set to a yoltage proportional to 
the pressure of the Siete gas in the cell with the switch 
closed to position (a). NH,D gas was obtained by mixing 


eeeeco® | ratio of “he and a commercially prepared sample of 


ND, in the cell. The gas was allowed to reach an equilibrium 


3 
condition before the relative intensity measurements were 
taken. Figure 26 shows the experimental relative intensity 


of the P(20) transition vs NH.D pressure. 


2 
An x - y plot of the relative modulated intensity 
vs the applied Stark cell bias voltage should provide the 
Stark shifted absorption lineshape. The experimental 
arrangement is the same as Figure 23 except that the D.C. 
Stark cell voltage is voltage-divided and fed to the x axis 
of the plotter with the switch in position {(b). The laser 
is set to oscillate on the P(20) line and monitored. The 
NH.D in the Stark cell is mixed in a 1:1 ratio as described 
above. The applied voltage to the cell was manually scanned 
from 0 - 3000 volts. Figures 27 - 34 show the experimental 
relative intensity vs applied Stark cell voltage for several 


NH.D pressures and scan directions. An unexpected discontin- 


2 
uity evident on the graphs and while monitoring the pyroelectric 
detector output directly on an oscilloscope showed up at 
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pressures above 5 torr. This abrupt change in the modulation 
of the laser made frequency stabilization impossible. 
Repetition of the experiment resulted in similar curves 
with no significant differences. It is evident from the 
curves in Figures 29 - 34 that the scan direction signifi- 
cantly effected the shape of the curve and location of the 
discontinuity, that the discontinuity occurred at different 
applied voltages to the Stark cell for different NH.,D 
pressures in the cell and that a substantial amount of 
modulation of the beam (as much as 50%) could be achieved 
using NH, D and the Stark effect. 

Closer inspection of Figures 29 - 34 indicated that 
the curve began to repeat itself after a discontinuity 
occurred independent of the scan direction. Figure 35 is 
an experimental plot of the relative intensity of the P(20) 
laser transition vs voltage applied to the Stark cell age 
the curves for both scan directions (a — O —- 3000v and 
b — 3000 - Ov) on the same graph. It shows clearly the 
effect of scan direction on the plot and repetition of the 
curve shape after the discontinuity. 

This repetition effect was studied by reversing the 
scan direction of the voltage applied to the Stark cell 
after each discontinuity yielding the plots shown in Figures 
36 and 37. The curves did repeat themselves with basically 
the same shape and with the discontinuities occurring at 
different Stark cell voltages and relative intensities. The 


curves of Figures 36 and 37 could be retraced for both 
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directions of scan provided the voltages corresponding to 
the discontinuity were not reached. Thus the shape of the 
irradiance is independent of scan direction. 

In order to observe the operation of the dither 
mechanism an A.C. voltage was applied to one plate while the 
D.C. voltage on the cell was maintained at a voltage where 
sufficient modulation was observed but below the point 
where the discontinuity occurred. A 200HZ A.C. signal was 
supplied to one of the Stark plates from an audio oscillator 
amplified to 60 volts vms. The laser was set to oscillate 
on the P(20) transition and the chopper was removed from 
the experimental arrangement in Figure 25. The signal from 
the pyroelectric detector was fed directly into the 
oscilloscope. A 200HZ sinusoidal signal was observed. This 
200HZ sinusoidal signal was lost when the laser was tuned 
off the P(20) line by changing the applied voltage to the 
piezoelectric driver in the laser. A scan through the 
laser's spectral signature revealed the expected results 


of modulation of the P(20) and P(14) laser lines by NH.D. 


2 
A 200Hz sinusoidal signal was evident on the oscilloscope 
for both these transitions with the amplitude of the signal 
for the P(14) line being about 50% of the amplitude of the 
Seonal for the P(20) line for the preset Stark cell voltage. 
Thus the Stark cell was actually modulating the laser beam 
and could produce a suitable signal for phase detection if 
the discontinuity could be corrected. 

Attempts to determine the cause of the discontinuity 
were unsuccessful. The voltage across the cell was equal 
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to the supply voltage (0 - 3000VDC) with no detectable 
voltage drop existing that might have caused a discontinuity. 
No current could be detected through the cell which would 
have indicated lonization of the gas or cell breakdown. [In 
order to determine if perhaps the beam was being steered 
off the detector by a refractive index gradient in the cell, 
rather than the gas modulating the beam, a detector was 
scanned across the output of the cell when the applied 
voltage on the cell was at a point where it was able to 
modulate the P(20) line. The pyroelectric detector was 
mounted on a moveable table and scanned across the Stark 
cell output signal giving a relative plot of intensity vs 
position. The signal from the detector was amplified, 
filtered and rectified and then fed into the y channel of 

a time base plotter. The scan rate of the detector was 
fixed so that the time base could be translated into a 
position measurement. Figure 38(a) shows the beam shape 
out of the cell as a function of position for zero voltage 
applied to the cell and Figure 38(b) the beam shape out of 
the cell when the beam was being modulated with an applied 
Stark cell voltage. The intensity of the beam is reduced 
with no noticable shift in position of the beam. The 
detector was positioned at the point of maximum intensity 
and a D.C. bias voltage was applied to modulate the base. 
The detector then moved to determine if the beam was being 
steered to the left or right, up or down. Figure 39 is a 


graph of intensity vs time and shows that the beam was not 
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EXPERIMENTAL PLOT TO DETERMINE IF A REFRACTIVE 
INDEX GRADIENT EXISTED IN THE STARK CELL. 


RELATIVE INTENSITY 





POINT ON GRAPH EXPLANATION 


Stark Cell Voltage Equals Zero 


tark Cell Coltage Increased to 
Modulate Laser Beam (Intensity Decreased) 


G Detector Moved to Determine if Beam 
was being Steered 

d Stark Cell Voltage Increased to Point of 
Discontinuity 

e Discontinuity (Intensity Jumps to a Zero 


Modulation Condition) 


FIGURE 39 
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being steered. The intensity as the detector was moved did 
not exceed the modulated intensity. Again the discontinuity 
is evident as the applied voltage was increased. 

The movable detector was placed at various positions 
in Figure 22 to determine if the beam shape was due to 
reflections in the cell or prior to the cell itself. Figure 40 
shows the relative intensity beam profile of the laser at 
three positions prior to the cell: 40(a) directly out of 
the laser prior to beam splitter No. 1, 40(b) between beam 
splitters No. 1 and No. 2 and 40(c) between beam splitter 
No. 2 and the cell. It appears that the germanium beam 
splitters cause some distortion of the beam and contribute 
to the distorted shape of the beam out of the Stark cell. 

The profile to the left of the main profile in Figures 38 (a) 
and 38(b) was probably caused by reflection from back surface 
of the parallel sided Ge beam splitter since it appears to 


be a mirror image of the main profile. 
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3. New Stark Cell Design 


In order to determine if the discontinuity was due to 
the nonuniformity of Stark cell plate spacing in the cell or 
the cell design itself, another Stark cell patterned after 
one described in Reference 31 was built. The new Stark cell 
design with its specifications is shown in Figure 41. Its 
design was an attempt to reduce the applied voltage require- 
ments by reducing the plate spacing to 3mm and eliminating 
the need for spacers to keep the plates parallel by attaching 
the plates to plexiglas holders with screws. This design 
allows for removal of the plates for modification, cleaning or 
substitution of a new design for the plates. The plate thick- 
ness can be used to vary their separation in the cell. The 
edges and ends of the plates were rounded to help prevent 
breakdown. The end caps hold the NaCl windows and simplify 
holding of the windows when modifying the cell. The cell held 
a good vacuum. The cell was filled with approximately 2 torr 
of NH,D with breakdown and ionization of the gas occurring at 
1000 volts. The breakdown was observed to occur along the 
edges and especially at the ends of the plates. It is believed 
that the inability to machine round edges is the primary cause 
for breakdown. Investigations into fabrication of plates with 
round edges, insulating the edges, use of geometrics other than 
strictly flat plates and electroplating of the plates on an 
insulator need to be conducted to produce a viable cell for 
frequency stabilization using Stark cell techniques. The pre- 


sent cell can be used to test plate designs and fabrication 
@echniques without having to build a new cell each time. 
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NEW STARK CELL DESIGN WITH INTERCHANGEABLE PLATES 


Nacl VOLTAGE PARALLEL GAS 
WINDOW CONNECTOR PLATES INLET 


SPECIFICATIONS 


PLATE: 

MATERTAL ALUMINUM 
LENGTH 15.0 cH 
WIDTH ee 
SEPARATION 3.0 mn 


CELL: 
MATERIAL PLEXIGLAS 





FIGURE 41 
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Iv. CONCLUSIONS 


A stable high voltage power supply (Sylvania 9415S) run 
from a stable line regulator was found suitable as a power 
source for the Sylvania CO, Model 941 Laser ina Stark cell 
stabilization system. 

A servo loop integrator was designed and built allowing 
the correction signal to the piezoelectric driver to be 
proportional to the time integral of the error signal pro- 
duced by phase sensitive detection. The integrator used in 
the feedback system should provide increased stability. 

A method for scanning through the laser's spectral signa- 
ture to find the required P(20) transition incorporating the 
Same integrator used in the feedback stabilization system 
was demonstrated and discussed. 

A Stark cell was tested to determine its characteristics 
for use in the stabilization system. It was found that a 
discontinuity existed in the cell's lineshape function which 
made it unsuitable for use in a phase detection stabilization 
scheme. The cause of the discontinuity could not be 
determined. A voltage drop across the cell, ionization of 
the NH5D gas, and beam steering caused by an index gradient 
being created in the cell were eliminated as possible causes. 

A second new Stark cell with the capability of modifying 
or changing the cell plates was designed and built:to deter- 


mine if the discontinuity found in the first cell was inherent 
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in the cell design or caused by a factor outside the cell. 
Breakdown occurred in the second cell prior to the field 
strength required for modulation making the second cell 
unsuitable for use in the system. The fabrication of the 
plates in the new cell was the major cause for breakdown. 
Further investigation needs to be conducted in the Stark 
cell plate design and fabrication to prevent breakdown 
inside the cell. 

Production of an error signal uSing a pyroelectric 
detector was demonstrated using the first Stark cell but 
was unsuitable for phase sensitive detection due to the 
discontinuity in the lineshape. It is expected that the 
successful fabrication and test of a Stark cell will produce 
a Stark stabilized co. laser in the assembled system. This 
system frequency stabilized on the P(20) transition is 
expected to be operational in both the laboratory and the 


field. 
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APPENDIX A 


MODEL 941 SYLVANIA LASER SPECIFICATIONS 


Output wavelength: 10.6 4m 

Power output: 3 watts 

Beam divergence: 4 mrad 

Beam diameter to 1/e point: 4 mm 

Beam polarization: vertical 

Cavity length: 45 cm 

Bore diameter: 6 mm 

Cooling requirements: tap water @ 15 to 30 ml/s 


mmeemc marror: 70% reflecting, dielectric flat 


Peezoelectric mirror: 100% reflecting, 300 cm radius of 
Curvature 
Piezoelectric -O4 parameter as given by 





QV 
Sylvania: -4 x ace mM vot 


Laser voltage: 14 kV @ 7.6 ma 

Tube voltage: 9.5 kV @ 7.6 ma 

Series ballast resistance: 600 kf 
Equivalent tube resistance @ 7.6 ma: 1.27 M 
Pressure broadened linewidth: 100 MHz 


Doppler broadened linewidth: 50 MHz 


al 
2L 


Cavity linewidth: 19 MHz 





mode separation: 333 MHz 
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APPENDIX B 
OPTICAL ENGINEERING —— LASER SPECTRUM ANALYZER 


The co. Laser Spectrum Analyzer is an infrared grating 
spectroscope designed simultaneously to monitor and visually 
display through the use of a calibrated thermal sensitive 
screen which darkens at points corresponding to the wave- 
length at which the laser is operating. The ovtical design 
of the instrument is baSically an Ebert Mount using a plane 
grating and a spherical mirror. Folding mirrors are used to 
reduce the length and the design modified so that the display 
screen can be located near the top of the instrument. 
Incorrect positioning of the instrument so that the laser 
light does not enter within the acceptance cone, will not 
create incorrect readings but simply cause the light to 
miss some element in the optical path and not reach the 
display screen. 

Bigure 1B is a side-view of the spectrum analyzer 
showing its essential components. In normal operation 
the portion of the laser beam which passes through the 
slit will strike a flat mirror which folds the optical 
beam and introduces a slight vertical angular component to 
the direction of the beam. The light then strikes a concave 
mirror which collimates the divergence introduced by the 


@iftfraction of the slit. The light then is reflected to 
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the diffraction grating. The grating is rigidly mounted 
on the flat mirror and held at a fixed angle. Only the 
light which is diffracted into the first order of the 
diffraction grating continues on to strike the curved 
mirror for a second time which focuses the image on the 
display screen. Each laser wavelength leaves the grating 
at a unique angle and is imaged at a different point on the 
calibrated scale. 

The specifications for the spectrum analyzer are given 


in Table I-B. 


OS 





Ind FEN 


*sgeT eTquqgentpy (01) ‘smazos yueuqsenlpe 
UTPEIQTTSO (6) ‘duet yeToTaezytn (g) feTv0E pexerqTTeD (L) ‘userog 


AvejTdsyq eaytzTeusg TeUIeyuL (9) 


‘BUTZBIF UOFZORIJITC (G) ‘torr TU 


eABDUND (7) ‘IOTITH 38TH (€) ‘tOzTUOW Meeg (Z) ‘ITIS STAPTIEA (1) 





ro La d+ 77 


6 8 


rezkTeuy wumr1z0edg JO META EPTS 


Or 





TABLE I-B 
SPECIFICATIONS 
The following is a list of the specifications for the 


Model l6A co" Laser Spectrum Analyzer. 


Wavelength range: 9.1 to 11.3 
Resolution: FOO see 3 ie) 

Maximum power: 75 watts 

Minimum power (focused): 0.01 watts 
Screen responce time: % second 


Grating efficiency: 90% horizontal polarization 
30% vertical polarization 


Optical design: 3/4 meter modified Ebert 
Total optical path length: 2.3 meters 
O O 
Acceptance angle: 2.5 horizontal, 1.5 vertical 
Slit width: adjustable - closed to 0.6 mm 
Wetont: 12 lbs. 
Mength: 18" 
PMiadth: 9.9" 
Minimum Mewes Ol. 3” 
Elevation of slit: 2.5" to 3.8" (adjustable) 
Electrical power: 115 V, 60 Hz, 10 watts 
2 , 
. c*°o, vibrational-rotational lines S O 
Ldentiried=s= 00 1 — O20 
P2 to P62 and R2 to R60 
00° = 10-0 
2 to P56 and R2 to R60 
Ol 


1 - 11°70 
P19 to P45 
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APPENDIX C 
SERVO LOOP INTEGRATOR 

The output of an ideal integrator is a voltage pro- 
portional to the integral of the input voltage; or in other 
words, the output is proportional to the product of the 
amplitude and duration of the input. The integrator performs 
this mathematical operation on an instantaneous basis. The 
result is an output exactly proportional to the area under 
a wave form. 

The circuit in Figure B-l performs integration using a 
MA 747 dual frequency compensated operational amplifier to 
force the same current through both RS and C._. The voltage 


if 
across the feedback capacitor is related to capacitor current by 


L fig dt (ke) 
S Cr 





< 
Hl 


Since the circuit causes i. to equal tne input current 
Vv. 
e — 1 


- f 
VN RB 8: ee Veit (2C) 
e Oo Ry Cer 1 


L 
The gain of the circuit is given by /RiCg- Thus the output 


a volts 
/C Ry 


The bias voltage, V5 and variable voltage offset null 


voltage will change by - /sec for each volt of output. 


resistor, R the drift inherent in the integrator. The 


V’ 


switch, Sy will reset the integrator to zero. 
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SERVO LOOP INTEGRATOR CIRCUIT 


SPECIFICATIONS ; 


OPERATIONAL AMPLIFIER: .4A 747 DUAL FREQUENCY COMPENSATED 


R, = R = Ry =f pa dee 
Ro smo k 0. 

Re = pelec@! 

Cn = Ty 

R= 10 (variable) 
Gain =~ /R Cp = = 

V5 = input voltage 

ae = output voltage 
Vi = bias voltage 

S. = reset switch 


Ky 


FIGURE C1 
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